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Abstract 

We give the description of discretized moduli spaces (d.m.s.) M'^g^n introduced in [Q| in terms 
of discrete de Rham cohomologies for moduli spaces A4g,n- The generating function for intersection 
indices (cohomological classes) of d.m.s. is found. Classes of highest degree coincide with the ones 
for the continuum moduli space Aig,n- To show it we use a matrix model technique. The Kontsevich 
matrix model is the generating function in the continuum case, and the matrix model with the 
potential iVatr {-^AXAX - \ log(l - X) - \X) is the one for d.m.s. In the latest case the effects 
of Deligne-Mumford reductions become relevant, and we use the stratification procedure in order to 
express integrals over open spaces M^g^^ in terms of intersection indices, which are to be calculated 
on compactified spaces ^/Cg"n ■ We find and solve constraint equations on partition function Z of our 
matrix model expressed in times for d.m.s.: = tr appears that Z depends only on 

even times and Z\t'^\ = C(qA^) e e ^*^'2„})+^({-*2"„>\ where is a logarithm of the partition 

function of the Kontsevich model, A being a quadratic differential operator in 
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Introduction. 



It has been shown recently that there exists a close but still not properly understood connection between 
three items: geometrical invariants of moduli spaces of algebraic curves; matrix models; integrable systems 
related to these models. The first relation was established by M.Kontsevich in Q who found a matrix 
model providing a generating function for intersection indices (integrals of the first Chern classes) on 
moduli spaces of algebraic curves. 

In this paper we describe some newly found applications of matrix models to the description of 
geometrical properties of the moduli spaces of algebraic curves. Here we should first mention brilliant 
papers by Maxim Kontsevich ||^ , in which the Kontsevich matrix model was introduced as the generating 
function for intersection indices of the first Chern classes on the moduli (orbi)spaces Mg.n of the surface 
of genus g and n punctures: 



Here uji is a Chern class associated with ith puncture. 

Kontsevich's papers were motivated by Witten's consideration |Q on two-dimensional (topological) 
gravity, or quantum gravity. Two different matrix model approaches were elaborated in order to describe 
such theory. 

The first approach is due to and it concerns a usual 1-matrix hermitian model with an arbitrary 
potential. In a "fat graph" technique starting with each graph we can construct the dual one corresponding 
to some Riemann surface with singularities of curvature concentrated in vertices of this dual graph. Then 
faces of this graph correspond to vertices in the matrix model graph and vice versa. If the initial potential 
contains only three valent vertices we can speak about "triangulation" of the Riemann surface. In what 
follows we shall deal with potentials of an arbitrary order, but we use the same term "triangulation" . The 
model with an arbitrary potential was solved exactly in |^ in the double scaling limit when the number 
of triangles tends to infinity and these singular metrics approximate "random metrics" on the surface. 
This model was presented by a hermitian N x N one-matrix model 



where P{X) = being times for the one-matrix model. For such system discrete Toda 

chain equations holds with an additional Virasoro symmetry imposed |5|. In the double scaling limit 
(d.s.l.) N oo and P{X) transforms in a way to incorporate surfaces with infinitely growing number 
of partitions and, as a result, the Korteveg-de-Vries equation arises. The partition function of the two- 
dimensional gravity for this approach is a series in an infinite number of variables and coincides with the 
logarithm of some r-function for KdV hierarchy. 

The second approach is based on cohomological considerations. In two-dimensional quantum gravity 
we have to integrate over space of riemannian metrics on manifolds modulo diffeomorphisms. Therefore, 
a finite-dimensional moduli space of conformally nonequivalent metrics arises. Integrals over such spaces 
have a cohomological description as an intersection theory on the compactified moduli space of complex 
curves. 

A fat graph technique was used in order to introduce coordinates on the moduli spaces. The coordi- 
natization means that we assign lengths U to all edges of the fat graph and the number of punctures, n, 
is the number of faces of the graph. We call this space Aig°^''. 

The model proposed by Kontsevich is a generating function for intersection indices or integrals of first 
Chern classes on the corresponding moduli space. It was a proposition by Witten ^ that these integrals 
yield correlation functions for the two-dimensional gravity coupled to the matter. In the continuum case 
the relation ( |3.l| ) holds where the integral goes over properly compactified moduli space and uJi are closed 
two- forms that are representatives of the first Chern classes of line bmidlcs on Mg^n- 




(0.1) 




(0.2) 
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For a general oriented graph of genus g and number of faces n the total number of edges (for trivalent 
vertices of the general position) is 617 — 6 + 3n which exceeds the dimension of Mg,n by n. So there are 
n extra parameters which are not related to the coordinates on the original moduli space itself. Namely, 
they arc perimeters of the faces of the graph. In the continuum case, due to Strebel theorem we 
have an isomorphism M.g^n ® R-" — ■^g"^'' and we define a projection tt : Aig°™^ R" to the space of 
perimeters. The fibers ^^~^(p^,) of the inverse map are isomorphic to the initial moduli space Mg^n f^nd 
hence they all are isomorphic to each other. 

Intersection indices for continuum case are expressed via the Kontsevich integral / DX exp{ tr ^KX^ + 
\X^^ with an external (Hermitian) matrix A. It satisfies equation of KdV hierarchy in times i„ = 
(2ri — 1)!! tr A"^""^. It is an asymptotic expansion of the string partition function 

T{t) = CXp ^(exp^i„a ) , (0.3) 

3=0 \ n I g 

and it is certainly a tau-function of the KdV hierarchy taken at a point of Grassmannian where it is 
invariant under the action of the set of the Virasoro constraints: CnT{t) = 0, n > — 1 0, ||T^. One 

might say that the Kontsevich model is used to triangulate moduli space, whereas the original models 
triangulate Riemann surfaces (see e.g. pl). 

In our recent papers we have proposed and developed an approach to the discretization of 

an arbitrary moduli space of algebraic curve. The connection of these spaces to a matrix model was 
established and also it was demonstrated explicitly that in the limit where discretization parameter 
becomes small this matrix model goes to the Kontsevich one . 

The discretization of A^^"™'' is rather simple - we assume that all lengths of edges are to be integer 
numbers (probably zeros). Fixing perimeters we always have a finite number of admissible sets of edges 
and a finite number of possible base diagrams for fixed g and n. Putting together all these possibilities 
we get the union of points of the discretized moduli space MJ^. 

We now are also able to define another projection tt : ■Mg'^^'^ ^+l^p G2Z4. '^here all perimeters are 
strictly positive integers with even total sum and consider its fibers 7r~^(p*). They are, generally speaking, 
finite sets of points belonging to the initial moduli space Mg^n- These sets are no more isomorphic to 
each other. Moreover, among these points there are always points which correspond exactly to singular 
surfaces ("infinity points"). We assume that the space of singular surfaces is dMg^n = Mg,n — Mg^n- In 
the usual Teuchmiiller picture all these points lie at the infinity, but in what follows we should include 
them explicitly into the game. We introduce an analogue of De Rham complex on these spaces using 
finite difference structures instead of differential ones. There are the spaces we call discretized moduli 
spaces (d.m.s.). Also instead of C/(l)-bundles for continuum case we shall consider "Zp-bundles" over 
these spaces. Thus we can define cohomological classes for d.m.s. as well: 



(0.4) 



There is a unique (up to isomorphisms) closed moduli space Mg,n = 7r-i(p,)[-^g?™''] and an infinite 
series of nonisomorphic 7r^^(p*)[Al^"'^], but for all of them the relation 

((T-di ■ • ■ Td„))g = (Tdi . . . rd„)g (0.5) 

holds true for higher-order integrals, '^di = d = Zg ~ i + n. We shall present matrix model arguments 
in favour of this statement, but note that, due to possible nonzero curvature of covering manifold, 
((rdj . . . Td^))g may be nonzero even for "Y^ di < d, in contrast to the continuum case. 



On the L.H.S. of ( p.5[ ) ((rdj . . . Td^))g can be presented in a form similar to (0.1) but with all quantities 
being related to the discretized moduli space (d.m.s.). 

One note about the structure of the simplicial complex described by the matrix model is in order. 
We denote this complex Tg^n, M.g^n itself is a quotient of it by some symmetry group of a finite order: 



3 



■Mg,n = Tg,n/^g- AH the integrals one can define on Tg^n instead of Mg^n, moreover, we have strong 
arguments in favour of the assumption that 7^_„ can be manifold (which is impossible in Teichmiiller 
picture). 

In the Kontsevich parameterization the evaluation of the integrals over 7r^^(p*) and 7r^^(p*) can 
be reduced to the calculation of integrals over volume forms on above mentioned finite coverings Tg^n- 
The discretization means in this language that we introduce an equidistant lattice on 7g,n and while 
calculating the volume we merely count a total number of sites in this lattice and divide it by some 
product of pf: . . .p°", where X)r=i a* = <^ = — 3 + n is the total dimension of A4g,n- Note also 
that all Tg n are compact spaces without boundaries. The sum over all points of a lattice is equivalent 
to the sum of unit cubes attached to each lattice point. Then we get the true volume of Tg,n only if all 
these points are nonsingular points, i.e., points of zero curvature. Since it seems not to be the case for 
every g and n (but holds for the case of torus with one puncture), some of indices {{di . . . dn))g may be 
nonzero for ^ di < d. 

In order to find a connection between moduli spaces Mg,n and d.m.s. we use a matrix model technique. 

Generalization of the Kontsevich model — so-called Generalized Kontsevich Model (GKM) is 
related to the two-dimensional Toda lattice hierarchy. It originated from the external field problem 
defined by the integral 

Z[A;N]= [ DX exp {N ti- {AX - Vo{X))} , (0.6) 



where Vo{X) — X^n^"^^'"'^" some potential, i„ are related to times of the hierarchy. This model is 
equivalent to the Kontsevich integral for Vo{X) ^ trX^. To solve the integral (p.6D one may use the 



Schwinger-Dyson equation technique 1 14 written in terms of eigenvalues of A. The Kontsevich model 
was solved in the genus expansion in the papers , ||l5|| for genus zero (planar diagrams) and in [p^ for 
higher genera. 

Another explicitly solvable model was introduced . The Lagrangian of this model has the following 
form: 



Z[A]^ J DXexp(^Ntri^~^AXAX + a[\og{l + X)-X]yj , A= diag(e' 



e-""). (0.7) 



This model may be readily reduced t o (|0.6|) w ith Va{X) = —X^/2 + alogX. It was solved in genus 
expansion in ||l7|| , |l8| . It appears (see ||19|, |20||) that it is in fact equivalent to the one-matrix hermitian 



model (0.2) with the general potential 

oo 

P(X) = ^r„trX", 
n=0 

where times are defined by the kind of Miwa transform (77 = A — aA~^): 

T„ = - tr??-" - — (5„2 for n > 1 and Tq = tr logT^-^ (0.8) 
n 2 

One of the motivations for considering such model was Penner's one-matrix model J DX exp tr (log(l— 
X) + X) whose asymptotic expansion gives the so-called "virtual Euler characteristics" of moduli spaces 
of punctured Riemann surfaces. These are positive rational numbers, which may be non-integer due to 
the orbifold structure of moduli spaces. 

It was demonstrated in that (3.7) is a model that describes in a natural way the intersection 



indices for the case of d.m.s. The only complication is that this model does not present generating 



function for the indices (3.4) straightforwardly because of the contribution from reductions. Indeed, any 
matrix model can deal with only open strata of a moduli space. It was not essential for the case of the 
Kontsevich model since there the integration went over cells of the highest dimension in the simplicial 
complex partition of the moduli space A4g_n- All singular points are simplices of lower dimensions in 
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Mg^n and give no contribution to the integral. But in the case of d.m.s. integrals over simplices of all 
dimensions are relevant due to the total discretization, so the integrals over reduced surfaces give nonzero 
contribution, which we should exclude in order to compare with the matrix model. The way to do it is 
to use a stratification procedure (2l|, which permits to express open moduli space A4g,n via Mg^n and 
moduli spaces of lower genera. 



In paper |g_2[ the explicit solution to the model (3/7), or, equivalently, to the general one-matrix model 



was found in genus expansion. The key role in this consideration was played by the so-called "momenta" 
of the potential resembling in many details "momenta" that appeared in the genus expansion solution 
to the Kontsevich model . We shall use some proper reexpansion of these momenta in terms of new 
times = tr c>^±i '^^ich stand just by the intersection indices (|0.4|). 

In the present paper we succeeded in finding and solving constraint equations for the model (l.C) in 
terms of times Tj^j. It appears that the model is readily expressed as a product of two Kontsevich models 
in times tj„^]^ = ±(2n + l)!!T'2^i intertwined by a mixing operator A that has the form of a canonical 
transformation expressed in terms of free fields of some conformal field theory. No limiting procedure is 
needed. 

The paper is organized as follows: 

Section 1 contains main notations and assertions of the paper. Also in Section 1 we solve constraint 
equations arising from Schwinger-Dyson equations in terms of times for d.m.s. We discuss the algebra of 
constraint equations and find the relation between the Kontsevich and the matrix model for d.m.s. For 
technical details of the proof, see Appendix A. In Appendix B one can find solutions of the constraint 
equations to few lowest orders of perturbation theory. A short review of the geometric approach to 
the Kontsevich model is given in Section 2. The definition of the discretized moduli spaces and the 
corresponding matrix model are presented in Section 3. Review of the previous results on explicit solutions 
of the Kontsevich, one-matrix model, and the model for d.m.s. is contained in Section 4. A detailed 
description of the simplest modular space, is contained in Section 5. Eventually, Section 6 contains 

a short summary of results and perspectives. 



1 Main results. 
1.1 Notations. 

Let g and n be integers satisfying the conditions 

5>0, n>0, 2~2g~n<0. 

Denote by Mg^n the moduli space of smooth complete complex curves C of genus g with n distinct 
marked points xi, . . . ,Xn and by A4g^n — a smooth compactification of Deligne-Mumford type. (The 
concrete scheme of this compactification will be discuss below.) 

Let Ci, i — 1, . . . ,n he line bundles on A^g,n. The fiber of Ci at (C, xi, . . . , Xn) is the cotangent space 
Let di, . . . ,dn be non- negative integers satisfying 

n 

di = diniQ ■Mg,n ^ ig — 3 + n, 

i=l 

and denote by (t^^ . . . Tci^)g the intersection index: 

^ n 

•'Ms.r. i=l 

where ci(£i) are first Chern classes of the corresponding line bundles taken in the moduli space A4g,n- 
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All matrix integrals are assumed to be integrals over Hermitian N x N matrices with the standard 
measure DX = Yli^j ddiXijd^Xij ni=i dXu. 

1.2 Main results. 

The main result by M.Kontsevich is the following 

Theorem 1.1. (M.Kontsevich j^) Considering matrix integrals over Hermitian matrices N x N as 
asymptotic expansions in times Tn = {n — l)!!tr pr+r we obtain 

CO oo ^ n 

g=0 "=i. s>o ^ ' ^di=3g-3+n i=l 

^ log ^r^- — = . . .), 

A = diag{Xi, Xn}. (1.1) 



It was E.Witten who first proposed that these intersection indices are nothing but correlators of 2D 
topological gravity. More, it was conjectured in and proved in |2^, |^ that these integrals satisfy 
equations of KdV hierarchy for times T„ . 

Consequently, in paper Q], we introduced a discretization of moduli spaces Alg,n and corresponding 
intersection indices 



J A4 1 



which can be non-zero even for di<d = 3g — 3 + n due to the discrete nature of this integral. 
One more complication is due to the necessity to integrate over the proper closure of the moduli space. 
In the continuum case it does not lead to any trouble since we may not concretize a conipactification 
procedure (one can choose the one due to Deligne and Mumford |2l|). The situation is different in the 
discrete case where singular curves give non-zero contributions to the integrals. Since there are no matrix 
model graphs corresponding to these singular curves we have to eliminate them using the stratification 
procedure which permits to "imitate" integrals over "open" part of discrete moduli spaces. 

Note, however, that apart from three-valent graphs, as in the Kontsevich case, we take into account 
explicitly graphs with vertices of arbitrary valence since they all contribute in the discrete case. These 
curves are non-singular. In this paper we prove the following 

Theorem 1.2. In the compactification scheme consistent with the Strebel parametrization of the moduli 
space Aig.n the following asymptotic expansion in terms of times 



is valid: 



oo oo ^ / j^xkil ? C I 

J2 (aiv)29-2+n E cg,„,,f--j 2^ "HI ^ryx 




rij rij 



k=l 




J £)^g-aWtr(iAXAX+ilog(l-X) + ix) ^ 
J jjX^-cNtriiAXAX^iX^) ^ •^A-p({T± }), 
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A = dmg{e^S.. ., e-^"}. (1.2) 

Here the sum runs over all reductions, Cg^n,rg cii^e positive rational numbers - coejjicients of reductions, 
\rq\ is non-negative integer, the power of reduction, < \rq\ < 3g — 3 + n, and, eventually, q is the number 
of components of the singular curve, q varies from 1 to 2g — 2 + n. Insertions of kj additional tq in the 
correlation function are due to reductions. 

The first term in (1.2), jr^j — 0, Cg_„^r, — 1 and q — 1 corresponds to the "highest" non-reduced 
modular space A^g*^^ from which we subtract integrals over reductions of the first power, which are again 
some closed moduli spaces of lower overall dimension, from which we are to subtract integrals obtained 
in points of intersections of these singular curves, et cetera. 

Our picture differs at this point from the one by the Deligne-Mumford compactification, where singular 
curve subspaces have symmetry groups of infinite orders. Therefore, in the Deligne-Mumford closure of 
the moduli space, these curves give no contribution to (1.2) since reduction coefficients Cg^n,rg are inversely 
proportional to the orders of the corresponding symmetry groups. It is not the case for the model (1.2) 
corresponding to the Penner-Strebel coordinatization picture. 

The matrix model on the R.H.S. of (1.2) was introduced in [|l^ where it was shown that it is exactly 
solvable in 1/iV expansion. It turns out p^ , [ pO| that with changed normalization factor it is equivalent 
to the Hermitian one-matrix model with general potential: 

J DXexp{~aNtT (iAXAX - i log(l - AT) - f )} = 

= (detA)-^-"^/2e-^''^^'e^''^(^+^")' / We-^(^), 



cx> ^ ^ 

C/(F) = 5:e.try", e» = ^tr (A + A-i)n -^^".^' (1-^) 

n— 1 ^ ^ 

where integral on the R.H.S. is done over Hermitian matrices of modified size ^ x ^ and we had to 
change the sign standing by the logarithmic term in order to keep this dimension positive. (In asymptotic 
expansion in 1/A^ one may easily make a transition back to "negative" dimensions by replacing N ^N.) 

As we deal with the general one-matrix model, we know that it obeys the equations of discrete Toda 
chain hierarchy in terms of times ^„ |^ . 

Note, however, that these times ^„ are of no use if we deal with times because their singularities 
are at different values of A. For ^„ they are A = i7r/2 -|- ink, k G Z, and for they are «7r(2fc -|- 1) for 
T+ and 2i7rfc for T", k E Z. 

Now we are going to discuss relations between the Kontsevich matrix model (1.1) and the model 
introduced in (1.2). We assume in what follows that the asymptotic expansion of (1.2) will be done in 
proper times T^^. An intermediate lemma states 

Lemma 1.2a. Partition function of the matrix model (1.2) J-kp{{T^}) depends only on even times 

-^2n — (2n+l)! dX^ c^±l 

The procedure of the double scaling limit (d.s.l.) permits to obtain the Kontsevich matrix model (1.1) 
starting from one-matrix model (l.A). Moreover, it is much easier to get (1.1) from the model(1.2) than 
from the one-matrix model. 

Namely, let us rescale integration variables in (1.2) as follows: 

a — > ae^^ 
X -*eX 

e^^e'^ (l.B) 
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In the limit e — > only those terms survive in the action of model (1.2) that give the Kontsevich action 
(1.1). It corresponds to the d.s.l. of one-matrix model with asymmetric potential. If we are going to 
consider the d.s.l. in the model with the symmetric potential where all odd are zero, then we have to 
choose a block-diagonal form of the external field matrix A = diag ( e '^^ , . . . , e ^"/^ , e "'"'^ , . . . , e ~-''«/2^ ^ 
Then in the limit e —^ we get two independent Kontsevich integrals over half-dimensional Hermitian 
matrices N/2 x N/2 taken with the same external field matrix A = diag (c^^^, . . . , e'*'"/=^). 

These relations were due to some limiting procedure which lead to some loss of information encoded 
in the model (1.2). It turns out that there exists a exact relation between models (1.1) and (1.2). In the 
present paper we prove the following theorem exactly solving the set of constraint equations: 

Theorem 1.3. Partition function of (1.2) and the Kontsevich model (1.1) satisfy an exact relation: 
where ^^„+i =±(2n+ iy.lT^„ and A is a quadratic differential operator in d/dT 



B2{n+m+l) 1 

^ 4(n + TO + 1) (2n + 1)! (2m + 1)! 

m,n— 



[J^ ^ , _A ^ _ 2f22("+"+i) - 1)^ — 



- 22»-i ^ 

',^^"(2^+1)! K^^'Wi^ 



C{aN) is a function depending only on aN that ensures that ^Kpi{T2n}) — 0, where = 0. i?2fc are 
Bernoulli numbers. 



1.3 Constraint equations for the model (1.2) 

In this subsection we present a sketch of the proof for Theorem 1.3 together with the algebra of constraints 
to the matrix model (1.2). (The complete proof is contained in Appendix A.) 

In what follows we treat all times (1.2a) as independent variables. Just like in the Kontsevich model, 
times {Tj^} (respectively, {TjT} become independent as iV — s- oo. However, there are mixing relations 
for these two sets that are valid (at least, formally) for all N. For instance, taking into account the pole 
structure of times (1.2a) in A variables we have 



tr ■ 



1 



= - tr 



1 



OO 



— = -}_^(fc + i)(z^)'^r,7 



fc=0 



These expressions always include infinite sums and, since the answer in (1.2) for finite g and n has finite 
polynomial structure, it is unique in terms of times {T^}. Therefore, in asymptotic expansion over TV 
and a we can treat all these times as independent variables. 



We begin with the matrix integral 



w{e^) = log 



'/ DX exp-aiVtr (iAATAX + i log(l - X) + iX) 
/ DX exp-aTVtr {\AXAX 



A 



(l.C) 



Integrating out angular variables we remain with the integral over eigenvalues Xi of the matrix X , for 
which we can write down the Schwinger-Dyson equations in terms of ( A.13| ). After a subtle algebra 



these equations can be reformulating in terms of times {T^f^} alone. The constraints acquire the form: 



k=0 



k-1 ' 



,u;(A) 



}-o. 



(l.D) 



k=0 
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Here 



22n+l 

(2n + 1)! 



2nB2n+i 



(— 

\2dX 



e-^± 1' 



where -B„(x) = X]r=o (n) ^s-^" ^^"^ Bernoulli polynomials, Bs being Bernoulli numbers. As we treat 
all times as independent, we get from (l.D) two independent sets of constraints on e™'-'*'^: 

2±e«'(A)^0, k>-l. 



Here L^^ is given by ( |A.36| ) and L±, s > 0, by ( |A.34| ). 

The constraint operators satisfy two halves of Virasoro algebra: 

[Lf,Lf]^^is-t)Lf^,, s,t>-l, 
[Lt, L-] - 0. 

Let us introduce creation-annihilation operators as follows (omitting ± signs): 

1 d 



2 dT- 



m > 0, 



m+7 



2m 
1 



1 + 77 T2,n, TO > 0, 



with corresponding commutation relations (for half- integer fJ.,iy E ^ + Z): 
The vacuum |0) is annihilated by with s < 0. 



We find the operator A (A. 39) of canonical transformation: 



— A + A 

e -^a„e-^ 
e-'^Z±e-^, s>-l, 



which completely split the dependence on "left" and "right" times in the Virasoro generators 

[aNf 



E 

m— — oo 



0^ 1 a* 



— -(aiV)a_,_^ 



+ 



(i.i?) 



Here the normal ordering : • : is defined w.r.t. the vacuum |0). 

These generators are nothing but Virasoro generators in the Kontsevich model (1.1). Therefore we 
get the assertion of the Theorem 1.3. 

One can interpret e'^^''^^2n+i})+-^A-({C2„+i}) jg) as a conformal vacuum | conf) of some c = 1 theory 
since it satisfies Virasoro conditions of the form: 



/:±|conf)=0, s>-l, 
where CJ^ obey two independent Virasoro algebra relations. 



(1.^^) 



A"] = 0. 



± are now related to the left and right moving sectors, which split completely on the level of Virasoro 
algebra. The Kontsevich integral makes a transition from a constant vacuum field |0) , which is annihilated 
by with ^ < 0, (/3s |0) = jg,^sfi for s > 0), to a conformal vacuum | conf), which satisfies left and right 
Virasoro conditions (l-F). 
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The operator A can be presented in an integral form. We assume that the numbers ±{aN) are related 
to eigenvalues of the momentum operators: 



P±\ 



Let us introduce a two-component bosonic field $(T, A) 



where 



.^(r.±, A) = ^ r±A"+i + a;± + p± log A + ^ 



A 



-n-l 



d 



n=0 



n=0 



n+1 dT^ 



± ■ 



Here the sum runs over all times, not necessarily even, derivatives act on the right. The corresponding 
currents, d(f){T'^,\) have normal ordering relations {d(f){X)d(l){iJ,)) ~ p^i^^. 

Then the operator A has the form: 



dX 



dfj, 



27ri / 2TTi 



where A(A + fi) is the following 2x2 matrix: 



A= (h — (h -^$^(T,A)A(A + /z)$(T,At), 



A(2/) = 



log 



l-e' 



+ isinh2y- 
'log(l+ e-J') 



log 



l-e 



log{l+e-y) 
- + isinh22/-|2/3 



(l.G) 



{l.H) 



This expression contains ambiguities. In fact, only symmetrical with respect to the change of variables 
y —y part of the matrix A(y) is rigidly fixed. The antisymmetrical part mixes odd and even time 



derivatives, and, therefore, gives zero when it acts on e 



^if({^2„})| 



The only non-zero contribution 



arises when mixing of occurs. (Two terms in the diagonal part, ^ sinh2t/ — and linear 



in y part of log ■ 



are combined in such a way that linear in derivatives term of (1.3) appears. 



2 Continuum Moduli Space Aig^n- 

2.1 Moduli space of algebraic curves and its pEirametrization in terms of 
Strebel-Jenkins differentials. 

We begin with an exphcit coordinatization of the moduli space Mg,n using the results of K.Strebel. He 
established the equivalence between "decorated" moduli spaces of algebraic curves and moduli spaces of 
ribbon ("fat", "oriented") graphs (R.Penner, J.Harer, D.Mumford and W.Thurston). 

A quadratic difFerential on a Riemann surface C is a holomorphic section of the line bundle (T*)*^^. 
In local coordinates it defines a flat metric on the complement of the discrete set of its zeros and poles: 

\^{z)\ ■ \dz\'^, where (p = ip{z)dz^. (2.1) 

All poles of the Strebel differentials are double poles in points of punctures. More, the quadratic 
residues in double poles are strictly positive real numbers. Since by the Riemann theorem for quadratic 
differentials # zeros - # poles = A{g — 1) (each is counted with its order), then for a general position 
point of Mg^n for the surface with n punctures there are 4(? — 4 -|- 2n simple zeros. 
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la lb Ic 

Fig.l. Horizontal and vertical lines of Strcbel metric in poles and zeros. 

A horizontal trajectory (geodesic) of a quadratic differential is a curve, along which ip{z)dz^ is real 
and positive. A vertical trajectory (geodesic), to the contrary, is the one along which (f{z)dz^ is real 
and negative. Let us consider the system of horizontal and vertical lines in the vicinity of the double 

2 

pole, where ipiz) = (2;r)^('^-z )^ i^^S- 1^)- It is easy to see that horizontal lines arc concentric circles 
around the point Zi while vertical geodesies are half-lines with the summit at the same point. If we take 
a fcth-order zero of f{z), then the situation is the following: there are exactly k + 2 horizontal and k + 2 
vertical half- lines with the endpoint zj (Fig. lb for fc = 1 and Fig. Ic for k = 2). For the general choice 
of the differential ip, horizontal trajectories are not closed, but Jcnkins-Strcbcl difFcrcntials are those 
for which the union of non-closed horizontal trajectories has zero measure. These trajectories are just 
those that connect zeros of the differential. Moreover, they decompose the surface into simply connected 
pieces (faces) with exactly one puncture in each. And, eventually, the lengths of all horizontal trajectories 
belonging to one face are the same - they are equal pi. K.Strebel proved the following: 

Theorem 2.1. For any connected Riemann surface C and n distinct points x\,...,Xn € C, n > 
0, n > ^(C) and n positive real numbers there exists a unique JS quadratic differential on 

C\{xi, . . . , Xn} whose maximal ring domains are n punctured discs Di, xi G Di, with circumference pi. 

The collection of these non-closed horizontal lines is an oriented graph corresponding to the Riemann 
surface. Orientation means that one may think about the edges of these graphs as strips with two sides 
belonging to two faces of the surface separated by this edge. (Looking forward, in matrix model technique, 
the index i will run along boundary line of a face.) 

Thus, in the Strebel metric each face converts into half infinite cylinder with the boundary consisting 
of borders of strips of the fat graph, i.e., the boundary is a polygon with the perimeter pi for ith face. 

What are the coordinates on the moduli space Mg,n in this picture? Let us supply the fat graph with 
additional data: we assign to each edge a positive real number Ig e R+, s being the number of the edge. 

These Ig have the meaning of lengths of edges of the fat graph T^, of a genus g. Is define coordinates on 
Mg"^'' - Qg — Q + 3n-dimensional linear space of graphs. Then pi are sums of lengths of edges incident 
to ith cycle. 

Zeros of JS differential correspond to the vertices of the graph, the valence of the vertex being equal 
to fc + 2 for a fcth-order zero. For a general case when all vertices are three valent the number of edges is 
equal to 6.g — 6 + 3n, which is the real dimension of the moduli space Mg^n plus n additional parameters - 
the perimeters of the faces. In order to find coordinates on Mg,n itself we must get rid of the dependence 
of these perimeters. 

Now we consider a set of all graphs with fixed g and n endowed with metric described above. To 

any JS differential we associate some graph. The inverse statement is also valid: having a graph one 
can construct the Riemann surface endowed with JS differential structure whose residues arc squares 
of perimeters of cycles on the graph. The set of all graphs modulo symmetry groups of graphs is a 
combinatorial moduli space A4g"J^^- The following statement holds true: 

Theorem 2.2. Let M^^^ denote the set of equivalence classes of connected fat graphs with metric 
and with valency of each vertex greater or equal 3. The map Aig^n ^ R" which associates 

to the surface C and positive numbers pi, ... ,Pn the critical graph of the canonical JS-differential, is 
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one-to-one. 

Thus considering all graphs of the fixed g and n we obtain a stratification on Ai^"™^ with the dimen- 
sions of strata equal to the numbers of edges. The open strata correspond to three- valent graphs and 
have the dimension 6g — 6 + 3n. 

We conclude this subsection with some notations from the graph theory necessary for what follows. 
For a fat graph F let Xq denote the set of vertices, Xi - the set of edges together with orientations 

defined for each edge, and X2 - the set of faces of the graph. Let sq and si be two permutations of Xi: 
Si changes the orientation of all edges simultaneously and is an involution, = id. Sq is defined as a 
rotation (clockwise, due to orientation) of edges incident to a vertex. 

Note: These transformations sq and Si are generators of a cartographic group corresponding to a chosen 
graph. The complete cartographic group can be represented as follows: let all edges be divided into 
two halves, all these halves being numerated. Therefore, we have a finite set of 2 x 7^ edges elements. 
The transformations Si, i = 1,2 define a permutation group on this set. These permutations do not 
necessarily preserve the edges of the graph as a whole, thus the authomorphism group of the graph is 
always a subgroup of the cartographic group and only in very few cases these groups coincide. 



2.2 Geometry of fiber bundles on 7Wg,„ and matrix integral. 

M.Kontsevich proposed a procedure for finding intersection indices (or, equivalently, integrals of the first 
Chern classes) on the moduli spaces. 

Let us consider a set of line bundles d whose fiber at a point E G Mg.n is a cotangent space to 
the puncture point Xi on the surface S. The first Chern class ci(£i) of the line bundle Ci admits 
a representation in terms of lengths of the edges Ij. The perimeter of the boundary component is 

Pi — "^laeii 

The first step in constructing ci(£i) is to determine a^, which is Z7(l)-connection on the boimdary 
component corresponding to the ith puncture. Since we already have an explicit coordinatization of the 
moduli space, we need only to make a proper choice of this connection in terms of Ij . It is convenient 
to introduce "polygon bundles" for each face - BU{1)'^°J^^ in Kontsevich's notations. These polygon 
bundles are sets of equivalent classes of all sequences of positive real numbers h,. . . ,lk modulo cyclic 
permutations. 

BU{iy°"^^ is the moduli (orbi)space of numbered ribbon graphs with metric whose underlying graphs 
are homeomorphic to the circle. There is an S'^-bundle over this orbispace whose total space EU{iy°"^'' 
is an ordinary space. The fiber of the bundle over the equivalence class of sequences . . . , is a union 
of intervals of lengths h . . . ,lk with pairwisc glued ends, i.e. a polygon. The inverse images of S'-'^-bundles 
are naturally isomorphic to the circle bundles associated with the complex line bundles d. 

Let us now compute the first Chern class of the circle bundle on BU (i)com5^ rpj^g points of EU (1)'=°"*'' 
can be identified with pairs {p, S) where p is the perimeter and 5 is a nonempty finite subset (vertices) 
of the circle R/pZ. Let <^i, < < . . . < (f)k < p,he representatives of points of S. The lengths of the 
edges of the polygon are 

k = 4>i+x - (i = l,...,A;-l), Ik = P + (l>i - (l>k- (2.2) 
A convenient form for 5^-connections on these polygon bundles is provided by 1-form a on EU{lY°"^^: 

k 

P \P 



1=1 

comb , ViTT(^ \comb 



a is well-defined and the integral of it over each fiber of the universal bundle EU{iy""^" —^ BU{iy 
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is equal to —1. The differential da is the pullback of a 2-form to on the base BU{1) 



comb 



l<i<j<k-l 

Extrapolating these results to the compactified moduli spaces we obtain that the pullback uji of the form 
uj under the ith map Mg_n x R-+ BU{iy°"^^ represents the class Ci{£i). 

Let us denote by tt : Mg°™'' R" the projection to the space of perimeters. Intersection indices are 
given by the formula: 

(r<i....r,J= / Yluf^, (2.5) 

J^-Hp^) t=i 

where = (pi, . . . ,p„) is an arbitrary sequence of positive real numbers and 7r~^(p*) is a fiber of Aig^n 
in M"""^^ 

We denote by the two-form on open strata of A4g°™^: 

n 

f} = ^p?^i(A), (2.6) 

1=1 

whose restriction to the fibers of tt has constant coefficients in coordinates (/(e)) . Denote by d the complex 
dimension of Aig^n, d ~ 3g — 3 + n. The volume of the fiber of tt with respect to fl is 



n 



d 



vol(7r-i(pi...,p„)) = / -W = T\ {plci{Ci) + ...+plci{Cn)) 

= E Xl^-JT^ra,...r,^)g. (2.7) 

One important note is in order. It is a theorem by Kontsevich that these integrations extend continu- 
ously to the closure of the moduli space Mg,n, following the procedure by Deligne and Mumford |2^]. (It 
means that we deal with a stable cohomological class of curves.) It is obligatory to consider a closure of 
the moduli space because all intersections can be consistently defined only on compact spaces. But it does 
not change insomuch our consideration until we integrate over continuum moduli space where nonzero 
contributions are given only by integrations over higher dimensional cells. All additional simplices which 
we add in order to close Mg,n are of lower dimensions. However, they will play a crucial role in the case 
of discrete moduli spaces. 

In order to compare with a matrix model we should take the Laplace transform over variables pi of 
volumes of fibers of tt: 

/•OO 

dp,e-P'^'pf' = (2d,)!A-2*-\ (2.8) 



for the quantities standing on the right-hand side of ( p.7| ). On the left-hand side we have 

/ .../ dp,A...Adpr,e-l^P'^' _ e", (2.9) 

Jo Jo JMg,„ 

and, due to cancellations of all pf multipliers with pi standing in denominators of the form fJ, we get: 

e^dpi A...Adpn = p Y[ die- (2.10) 

eeXi 

We use standard notations: Xq is the total number of g-dimensional cells of a simplicial complex. {Xi 
is the number of edges, Xq ~ the number of vertices, etc), p is a positive function defined on open cells: 

P^\I[\dp^\^^) ■■ n (2-11) 
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Surprisingly, the constant p in fact depends only on Euler characteristic of the graph F, p = 2 

^ ^ ^d+#x,-#Xo _ (2.12) 

The integral 

IgiK):^ f exp{~J2x,p,) l[\dl{e)\ (2.13) 

is equal to the sum of integrals over all open strata in Al™™''. These open strata are in one-to-one 
correspondence with a complete set of three-valent graphs contributing to this order in g and n. It is 
also necessary to take into account internal automorphisms of the graph (their number, in fact, shows 
how many replica of moduli space one may find in the cell) . The last step is to present the sum ^ \iPi 
in a form dependent on 1^: 

n 

^A,K= ^ /e(AW+A(2)) (2.14) 

1=1 eSXi 

for each graph. Here A^^^ and X^iP are variables of two cycles divided by eth edge. Performing now the 
Laplace transform we get the celebrated relation |^ (Theorem 1.1): 

di,...,<i„=0 i=l P TT V ; I J 

where the objects standing in angular brackets on the left-hand side are (rational) numbers describing 
intersection indices, and on the right-hand side the sum runs over all oriented connected trivalent "fat- 
graphs" r with n labeled boundary components, regardless of the genus, is the number of vertices 
of r, the product runs over all the edges in the graph and # Aut is the volume of discrete symmetry 
group of the graph T. 



The amazing result by Kontsevich is that the quantity on the right hand side of (2.15) is equal to the 
free energy in the following matrix model: 

, (A) „ / rfX exp aiV (- 1 tr AX^ + ^ tr X^) 

J dX expaN {^^ti-AX^) ' ^' ' 

where X is an A'^ x A'' hermitian matrix and A = diag (Ai, . . . , Aat). a is an additional parameter 
enumerating the boundary components. Though each selected diagram has quantities (A^ + Xj) in the 
denominator, when taking a sum over all diagrams of the same genus and the same number of boundary 
components all these factors are canceled with the ones from nominator. 

Feynman rules for the Kontsevich matrix model are the following: as in usual matrix models, we 
deal with the so-called "fat graphs" or "ribbon graphs" with propagators having two sides, each carries 
corresponding index. The Kontsevich model varies from the standard one-matrix hermitian model by 
additional variables Xi associated with index loops in the diagram, the propagator being equal to 2/ (A^ + 
Xj), where Xi and Xj are variables of two cycles (perhaps the same cycle) incident to two sides of the 
propagator. There are also trivalent vertices presenting the cell decomposition of the moduli space. 

Let us consider the simplest example of genus zero and three boundary components which we sym- 
bolically label Ai, A2 and A3. There are two kinds of diagrams giving contribution into this order (Fig. 2). 
The contribution to the free energy arising from this sum is 

1 If 1 



6(Ai + A2)(Ai + A3)(A2 + A3) 3 i 4Ai(A2 + Ai)(A3 + Ai) 

(1 ^ 2, 2 ^ 3, 3 ^ 1) + (1 ^ 3, 3 ^ 2, 2 ^ 1) 

2A1A2A3 + A2A3(A2 -l- A3) -l- AiA3(Ai + A3) + AiA2(Ai + A2) 
12AiA2A3(Ai + A2)(Ai + A3)(A2 + A3) 

I2A1A2A3 



14 



This example demonstrates the cancellations of (A^ + Aj)-terms in the denominator. 



1/6 




perm. 



Fig. 2. The g=0, s=3 contribution to the Kontsevich model 



The quantity standing in the R.H.S. of ( 2.15| ) is nothing but a term from 1/A^ expansion of the 
Kontsevich matrix model. Eventually, we have: 



2-2g^2-2g-n 



E 



3 = 
n=l 



si+2s2 + ...+fesfc=d 



'sol. ..Ski 



n 



tr 



(2d,-l)!! 

A2<i. + 1 



log- 



(2.18) 



Thus, the Kontsevich matrix model provides a generating function for the intersection indices of the first 
Chern classes on the moduli (orbi)spaces. 

Let us introduce a new important object - the 2-vector (3 which is defined on the higher dimension 
cells of A^g „: 

2 ~ 



d 



dl{x) dl{so{x))' 



(2.19) 



Here sg is an authomorphism of A^™™'' which "rotates" each edge x by |7r clockwise over a vertex. This 
2-vector defines a Poisson structure on the cells of the higher dimension in A/t™™''. In order to see it let 
us calculate its kernel: 



Proposition 2.1 Kcr/3 = 7r*r*R'J:. 

Proof. Ker/3 is a space of functions on the edges of the graph such that the following relation holds: 

h + h^h + h. (2.20) 

If we take all four edges neighbour to the edge /o, the order is clockwise. (If we combine all terms from 
(2.19) with d^ = we just get A (9i +8^,-82 — 84). Hereafter we shall denote the derivative over li 



by di. There exists a unique function g defined on the set of faces (boundary components) X2 such that 
on each edge fi — g'"^'^ + gf'^ - the sum of g-variables of faces incident to this edge. In the neighborhood 



of each vertex g can be reconstructed inambiguously: gi — (/2 + /s — /i)/2. Condition ( 2.2C ) ensures 
that moving from vertex to vertex this number is preserved. This f/-function does not coincide with 
perimeters pi for each face, but one may find the relation between them. (In order to construct pi from 
the set of g-variables one may simply take the sum of fi — + 5-^^ over all edges surrounding the face. 

Thus M.g_n is a Poisson manifold whose symplectic leaves are fibers of the projection tt. The following 
proposition again by M. Kontsevich establishes the relation between (3 and Chern classes tOi ( ^.4[ ): 



Proposition 2.2. Onn^^ijp^) AP^'^ = Pi 



X W,;. 



i=l 
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3 Discrete Moduli Space Mg^^ • 

We describe a discretization of the moduli spaces in this section. We hope that these discrete moduh 
spaces would be helpful when taking a quantum deformation of the Poisson structures on the moduli 
spaces. Also, these discrete spaces have their proper meaning because, as we shall show, they admit a 
nice description in terms of another explicitly solvable matrix model. We shall start with a description 
of this model, merely for being acquainted with it. 



3.1 The matrix model for d.m.s. 

Let us consider the matrix model (l.C) |l7j, where we denote, for simplicity, /ii = e"**'. It includes, 
in contrast to the Kontsevich model, all powers of X" in the potential since it describes the partition of 
moduli space into cells of a simplicial complex, the sum runs over all simplices with different dimensions. 
(In the language of the Kontsevich model the lower the dimension is, the more and more edges of the fat 
graph are reduced). 

The logarithmic potential makes this model similar to the Penner matrix model / DX exp tr (log(l- 
X) + X^ counting virtual Euler characteristics of the moduli spaces. 

We find the Feynman rules for the theory (l.C). First, as in the standard Penner model, we have 
vertices of all orders in X. Due to rotational symmetry, the factor 1/n standing with each X" cancels, 
and only symmetrical factor l/#Autr survives. Also there is a factor {a/2) standing with each vertex. 
As in the Kontsevich model, there are variables fii associated with each cycle. But the form of propagator 
changes — instead of 2/{Xi + Xj), we have 2/a{fiifij — 1). 

Let us consider the same case (g = 0, n = 3) as for Kontsevich model. One additional diagram 
resulting from vertex X^ arises (Fig. 3). 



a-V3 






-perm. 



Fig. 3. g=0, s=:3 contribution to the model (l.C). 



Symmctrization over fii, fi2 and /is gives: 

3 \ 2(/ii/i2 - l)(AiiAi3 - 1) 
1 r 2a-iiV2 



perm. 



2a-^N^ 



6(miM2 - l)(AiiA«3 - 1)(M2A'3 - 1) 



3 \ 2{^il-l){n^,l2-l){^^l^ii~l) 
Again, collecting all terms we get: 



perm. 



(3.1) 



2a A^ I V- _ 9 I M2A^3 - 1 - 1 MiA*3 - 1 ~ 

and after a little algebra we obtain the answer: 



2 -1 A*iAi2 + AiiAi3 + Ai2Ai3 + 1 



3(M?"l)(Mi-l)(A^i-l)' 



(3.2) 



(3.3) 
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We see that here, just as in the standard Kontsevich model, the cancellation of intertwining terms in the 
denominator occurs that leads to the factorization of the answer over l/{pf — l)-terms. This simplest 
example shows that there should be some underlying geometric structures in this case as well. 

Note that technically the reason why (l.C) depends only on tr A'^ {k < 0) is the following: This model, 
as well as the Kontsevich one, belongs to the class of Generalized Kontsevich Models jl^. It means that 
after some simple transformation we get from (l.C) the model with the potential AX + V{X), which 
depends only on Miwa's times. 



3.2 Moduli Spaces and discrete De Rham cohomologies. 

Let us now consider a discretization of the moduli spaces M.g,n and A^™™^. We shall consider the 
following set of parameters U : 

n 

k e z+ u {0}, e z+, e 2Z+. (3.4) 

1=1 

So all li and pi are now integers, but some of U can be zeros while all perimeters are strictly positive. The 
sum of all perimeters is even because each edge contributes twice into it. We call this (combinatorial) 
space A^*^*^. It is worth mentioning that now we explicitly include into play such points of the original 
Mg,n which are points of reductions (singular curves) and curvature (orbifold points that are stable under 
the action of some non-unit subgroup of the symmetry group in the Teichmiillcr space). While keeping 
all Pi fixed in a general case we can put a number of Ij exactly equal zero. Some of these configurations 
belong to the interior of A^g.n, but not all — it means that among the points of A^*^'^ there are points 
that lie on the boundary dAig^n. Such points correspond to the reductions of the algebraic curve. Also 
we shall use the notation ■Mg^n'^ for such subset of M.'^g^^ where all points of reduction are excluded. 



It turns out that this choice for d.m.s. is rather natural since all the quantities (2.3-2.6|) have corre- 
sponding contemplates in this discrete case. 

First we define the action of the external derivative d and the integration over these (orbi)spaces. We 
shall write the d-action on functions, the extrapolation to the space of skew symmetric forms is obvious: 

k 

df{h, . . . , /fc) = ^(/(?i, ...,k + l....M)- fih, h))dk (3.5) 

i=l 

As for the integral over domain fi, there is again a proper generalization of it to this discrete case: 

f{h,...M)dh...dh:= V fih,..., Ik). (3.6) 

Instead of BUiVf"'^^ we have (orbi)space of equivalence classes of all sequences of non-negative 
integers li,.. .,lk modulo cyclic permutations. An analog of S'^-bundle is now a kind of Zp~ "bundle" 
over this new discrete orbispace whose total space iJZ™™'' is an ordinary rectangular lattice. The fiber 
of the bundle over the equivalence class of sequences li, ... ,1k is again the polygon with integer lengths 
of edges h. . . ,lk. 

Let ipi be coordinates on -BZp, just as in ( ^.2[ ): 

= - 0i (i = 1, . . . , fc - 1), Ik = p + 4>i - 4>k. (3.7) 

Due to the linearity of (|]^) in li and (jjj it can be straightforwardly generalized to our case. Denote by 
a the 1-form on EZ^"™^, which is equal to 

a = ykx^. (3.8) 
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The integral of a over each fiber of the universal bundle ET^l"""^ is equal to -1. The 

differential da is the puUback of a 2-form uj on the base i^Z™'"^, 

LJ = > — A^. (3.9) 

l<i<j<k-l ^ ^ 

Extrapolating these results to the whole discrete moduli space we obtain that the puUback Ui of the form 
oj under the ith map Mg^ BZ™™'' represents the class ci{Ci). 

Let us denote by tt : ■M'^^'^ — > [Z"]^^^^ the projection to the space of perimeters with the restriction 
'^iPi G 2 • Z+. Intersection indices are given again by the formula: 

{{rd.---rdj,= L n^''' (3.10) 

where p^, — (pi, . . . ,p„) is an arbitrary sequence of positive integer numbers, J2^=iPi being necessarily 
even, and TT~^{p^) is an analogue of the fiber of A^g,n in -^g'rT- Note, however, that now the volume 
of Tr~^{p^) may depend on pi in a non-monomial way. Therefore, these indices may be non-zero for 
di < d = 3g — 3 + n. 

One important note is in order. Each fiber tt^^{p^,) contains finite number of points. Thus, these 
fibers are not isomorphic. But they all are analogues of the initial moduli space Mg.n taken with different 
perimeters. For this reason we call them discretized moduli spaces. It appears that relation (3.1C) remains 
valid independently of how many points from the initial Aig^n give contribution to the fiber tt^^{p*) (it 
can be even only one point of reduction, as we shall see for A^i.i). Values of these intersection indices 
are some rational numbers due to the orbifold nature of the initial moduli space Mg,n- This nature 
reveals itself as symmetries of the graphs. But these symmetry properties are the same, whatever case 
- continuum or discrete, and whatever values of perimeters we choose. Thus, preserving symmetry 
properties we preserve the values of cohomological classes on both continuum and discrete moduli spaces. 

The difference appears when we consider ((rdj . . . Td„ ))g with ^ < 3g — 3 + n. In the continuum case 
such quantities vanish in contrast to the discrete one where they may be non-zero due to the possible 
non-zero curvature of the covering manifold. 



3.3 Matrix Integral for Discretized Moduli Space. 

We denote by V, the two-form on A^g'^'^: 

n 

h = Y^p^iUi, (3.11) 

whose restriction to the fibers of n has constant coefficients in coordinates (^(e)). d is again the complex 
dimension of A4g,n, d = 3g — 3 + n. The volume of the fiber of tt with respect to is 
^ r fid If ^ ^ d 

V0l{n-^{pi . . . ,Pn)) = / -TT = 37 L {Pl^l + ■■■+ pi^n) = 

= E \{%{{^d,---^djg. (3.12) 



Next step is to do a Laplace transform in both sides of (3.12). Of course, now we should replace 
continuum Laplace transform by the discrete one and also explicitly take into account that the sum of 
all Pi is even. On the R.H.S. we have: 

li / o \ 2di if"! " 1 ~1 



^Pie2Z^ 



=1 
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On the L.H.S. of (3.12) we again substitute 



e"' 



dpi A . . . dpn\Y ^^^1 = p W die. (3.14) 



Here the constant p is the ratio of measures similar to (2.11) and we only need to take into account the 



restriction that the sum of all pi is even. It leads to the renormalization of p for the case of d.m.s.: 

p^p/2, (3.15) 



where p is given by (2.12) 



Now we give a matrix model description for these "new" intersection indices. Here we immediately 
encounter some troubles. Let us consider the correspondence between graphs and different points of 
7r~^(p*). First, there are points of a general position with all U greater than zero, which correspond to 
graphs with only trivalent vertices. Second, there are such points of TT~^{pi,) where some U are zeros, 
but these points still do not correspond to reductions. For example, see Fig. 4, where for the torus case 
Fig. 4a represents a point of a general position, > 0, i = 1,2,3, and Fig. 4b gives an example of the 
graph for which one (and only one) of U is zero. Such graphs do not correspond necessarily to singular 
curves. The one depicted in Fig. 4b determines the curve in M.i^i (in Teichmiiller parameterization) with 
purely imaginary modular parameter t. Of certain, if we want to include such graphs into consideration 
we should take not only trivalent vertices, but vertices of arbitrary order. In the continuum limit we 
did not take into account such graphs, since they correspond to subdomains of lower dimensions in the 
interior of the moduli space, the integration measure being continuous and we may neglect them. Here 
the situation is different and we should take into account all such diagrams as well. 





^2=0 





Fig.4a Fig.4b 

Fig. 4. Diagrams for different regions of A^i_i 



Fig.4c 



But in each 7f~^(p*) there always are (except for the case M.0,3) true points of reduction (see, for 
example. Fig. 4c when two of U are zeros). We are not able to give a matrix model description to such 
curves. At first sight it would mean that all the construction fails since we still do not discuss how to 



"exclude" such reduction points from tt ^(p*) by modifying somehow the relation ( 3.10 ). Let Aif'^'^ be 

such subset of A^*^'^ where all points of reduction are excluded. Thus we need to release the integration 
over open Mg,n from the total integration over Mg,n- In order to do it we use a stratification procedure 
a'la Deligne and Mumford [Q. The idea is to present the open moduli space Mg^n as a combination of 
Mg.n and Mgj^nj of lower dimensions. The description of this procedure one can find in | |Tl[ |. 

Geometrical meaning of the reduction procedure is that we subsequently pinch the handles of the 
surface (Fig. 5). One can see that there are two types of such reduction: for the first one, by pinching a 
handle we result in the surface of genus 5 — 1 and two additional punctures. Thus, from the space Aig,n 
we get after such reduction Mg-i^n+2 (Fig. 5a). In the second case by pinching an intermediate cylinder 
we get two surfaces of the same total genus and two new punctures: one per each new component. It 
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means that the initial moduh space Aig^n spHts into the product A^g^.m+i ^ -^92, "2+11 5i + .92 = <7, 
ni + 712 = n (Fig-5b). 




I 



Fig. 5a. One-component type of the reduction. 



pinching 




Fig. 5b. Two-component type of the reduction. 

One may easily check that total complex dimension of the resulting spaces is in both cases d — 1, 
where d — 3g — 3 + n is the dimension of Aig^n- So the general receipt of how to express A^g,n via closed 
moduli spaces is to construct an alternative sum over reductions: 

3(/— 3+n 

Mg,n^ E (-l)'^'®l=l>'9.,™.+fc., (3.16) 



where the sum runs over all q-component reductions, Vq is the reduction degree and kj being the number 
of the additional punctures due to reductions. The dimension of A4gj.nj+kj is dj = 3gj — 3 + nj + kj, 



q q 
^nj=n, ^dj^d-Tq. 



(3.17) 



Thus we have: 



dpi A ... A dp„ 



1 



dpi A ... A dpn + 



3f/— 3+n 



^ (-ir- ^ / Efa^^a ei:.^'^-dpiA...Adp„^. 



(3.18) 



educti 
'■0 = 1 



Now we can find, using (3.14), a matrix model description for the L.H.S. of (3.18). Just as in the 
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continuum case we have 
L 



H.S.= / expj- 5] ;e(AW+A(2))| xpx [] \dl{e)l 



(3.19) 



where p = 2''-+*^^^*^o-i ^ 

This last expression can be presented as a sum over all possible "fat graphs" V with vertices of all 
possible valences for given genus g and number of faces n. We should again take into account the volume 
of the authomorphism group for the graph. This volume coincides with the number of copies of equivalent 
domains of the moduli space Mg^m which constitute this cell of the combinatorial simplicial complex. 
Finally, we "integrate" over each Z(e), i.e. , take the sum over all positive integer values of l{e) (because 
we already took into account all zero values of l{e) doing the sum over all graphs). Eventually, we have: 



L.H.S. = 2"^-^ J2 



1 



#Aut (r) 



X 2" 



n 

eeXi 



Graphs r 

It is nothing but a term from the genus expansion of the matrix model (l.C) with 

A = diag{e^S . . . , e^"}. 
Then log Z [A] has the following genus expansion: 

oo oo 

logZ[A] = ^ ^(7Va)2-2<,^-„^_^(^^^ ^ ^^^^y 

3=0 n=l 



(3.20) 



(3.21) 



(3.22) 



Let us use the relations (|3.13| ) in order to express the R.H.S. of ( p.lSD via intersection indices. 

1 



< An) 



X tr 



n 



lk=l 

Theorem 1.2 is proved. 



d 
dXk 



2dt 



Crii 
t=i ""^"-^ 




(3.23) 



Formula (3.23) is our main result. For practical reasons it is sometimes convenient to rewrite 



E 




60 + ''l + --- + i'fc = '>3 
i)0 + l-''l + - - + '= i'fc=33j-3 + s 



(((ro)'"'...(r,)'"=ro...ro))3,. 



(3.24) 



Taking this expression for the case A^o,3 (without reductions) and recalling that (tq )o = 1 we immediately 
get the answer (5.3) after a substitution pi ~ e^K 

Since the matrix model (l.C) is equivalent to the hermitian one-matrix model with an arbitrary 
potential, formulae ( 3.22 )-( 3.24 ) above give the solution to such models in geometric invariants of the 
d.m.s. 



There are two complications in the final relation (3.23) that make it qualitatively different from the 
Kontsevich model. First of them is the "sum over reductions". Another is the "new" averaging ((. . 
The sum over reductions turns to be rather involved for the following reasons: When we considered the 
orbits 7r~^(p*) we assumed that they belonged to one copy of the moduh space M.g^n- But when we 
deal with the cell decomposition it is much more convenient first to consider the total simplicial complex 
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(which we shall denote and only afterwards take into account internal automorphisms of which 
eventually produce Mg,n as a coset over a symmetry group Gg,„. One may consider, instead of tt and 

TT, mappings /3 and ^, respectively, /3 : 7^,„ «) R'^ -> '^^^ ' '^g'^,n''\P*] ^ i'^+Uven ^ ^g^^: 

where T^^'^ are again finite (nonisomorphic) sets of points of 7g_„ supplied with the discrete de Rham 
cohomology structure. 

For these spaces 7^_„ an analogue of the formula ( |3.23 ) exists. The only difference is that we should 
multiply all indices {{rd-^ ...Td^))g by the order of the symmetry group Gg^„. This number of copies 
multiplied by the order of the group Gg-^nj+kj is not necessarily divisible by the order of Gg^n- Thus 
when we write in formula ( |3.23| ) "the sum over reductions" we bear in mind that the coefficients in this 
sum are not necessarily integers! (See Section 6 for an example). 

The next section is devoted to the comparison of the matrix integrals in the Kontsevich and the 
matrix model for d.m.s. using exclusively matrix model tools, which permits to prove the coincidence of 
(tjJj . . . Td^)g and ((rdi . . . Td^))g in the highest dimension X^ILi di — d = ig — i + n. 



4 Comparison of two matrix models 



This section is based on the results of papers |l8j and ||l2|. It was explicitly demonstrated in |19 
that the matrix model (l.C) is equivalent to the standard hermitian one-matrix model 



NxN 



d(t)exp{-NtrV {(!))), 



where the integration goes over hermitian NxN matrices and 

y(0)=f;^0^- 

is a general potential. Then the following relation holds: 

Z[g,N{a)] = e-^*"i'/^Zp[i^,N], N{a) = -aN. 
Here the partition function Zp[ri,N] is 



Zp [77, TV] 



dX exp 



NxN 



NtT ( -rjX- ^X^-a\ogX 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



the integral being done over hermitian matrices of another dimension NxN and the set of the coupling 

(4.5) 



constants (4.2) being related to the matrix rj by the Miwa transformation 



9k ^ ^^^'^ ^ for fc > 1, 50 = -^ tr \og-q ^ 



(Note the changing of sign in front of the logarithmic term.) Now after the substitution 

77= ^/^(A + A-l) (4.6) 

and after the change of variables X ^ {X — l)A^/a we reconstruct the integral (l.C) (with a multiplied 
by two). 

Note that we can do a limiting procedure (which is a sort of the double scaling limit for the standard 
model (4.1)) resulting in the Kontsevich integral (2.2) starting from the model (l.C). It looks even more 
natural in terms of this model than in terms of the one-matrix integral ( |4.l| ). Namely, let us take in 
(l.C) 

1 



A: 



(4.7) 



Then after rescaling X eX in the limit e ^ we explicitly reproduce (2.2) from (l.C). During this 
procedure we can keep the size N of matrices of (l.C) constant, but the size N{a) of the matrices of 
hermitian model goes to infinity in the limit e ^ 00. 
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4.1 Review of the solutions to Kontsevich model and model (l.C). 



Since the models (l.C) and (4.1) are equivalent, we can use the explicit answers for ( |4.l| ) found in |_2| 
in order to check the validity of our formulae (3.23) and to compare the values of intersection indices in 
both Kontsevich model (2.16) and the model (l.C). Both these models were solved in genus expansion 
in terms of the corresponding momenta. For the Kontsevich model this solution was presented in 
and for (14) or, equivalently, (11) — in ||l^. Here we present the results. (Throughout this section the 
expansion parameter a should be replaced by —2a in order to compare with the results of [p2|.) 

1. The solution to the Kontsevich model is 



l0gZK[iV,A]=^7V2-29^^ 
9=0 



Kont 

a 



For the genus expansion coefficients we have 



F, 



Kont 



E 



>l 

(Q.-l) = 3g-3 



IoL\ • • • 



for 5 > 1, 



(4.8) 



(4.9) 



fc > 0, 



where {■)g are just intersection indices and the moments /fc's depending on the external field M are 
defined by 

1 1 ^ 1 
~ {2k-iy.l N ^ (to2 _ 2uo)''+^/^ 

and uq{M) is determined from the equation 

Mo = Io{uq,M). 



(4.10) 



(4.11) 



2. The solution to the model (4.4) can be written as 

oo 

l0gZp[7V,77]=^iV2-23^^^ 

where 



(4.12) 



3=0 



F,= (ai...a.;/3,...Ak/3,7)3 ^"^--,^";/;:---'^'' 5 > 1- (4.13) 

This solution originated from the one-cut solution to the loop equations in the hermitian one-matrix 
model, X and y being endpoints of this cut, d — x ~ y, and for momenta A/fc, we have 

Mk =^E7—: I777-'5m fc > 0, (4.14) 



Jk 



{rjj — x)^+^/'^{rij — y)-'-/^ 



ly \ 



Jk,l 



fc > 0. 



(4.15) 



The brackets {■)g denote rational numbers, the sum is finite in each order in g, while the following 
restrictions are fulfilled: If we denote by Nm and Nj the total powers of M's and J's, respectively, i.e. 



Nm = s-a, Nj = l-f3, 
then it holds that Nm < 0, Nj < 0, and 

Fg : Nm + Nj^2~2g, 

s I 

Fg: 5](a.-l)+^(/3,-l)+7-45~4 
i=i j=i 

s I 

Fg: ^(a,-l)+^(/3,-l)+7<35-3 
i=i j=i 



(4.16) 



(4.17) 



23 



We again have nonlinear functional equations determining the positions of the endpoints x and y: 

1 ^ 1 
^ ~i ViVz - x){ri, - y) 



x + y 



0, 



1 ^ 
N ^ 



(a; - yf 



~t ViVi - - y) 



-2a + 1. 



(4.18) 
(4.19) 



The solutions to the first two genera have, as usual, some peculiarities. For 5 = 1 we have 

1 
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logMiJid*, 



(4.20) 



and for zero genus we have, after taking a double derivative in a in order to exclude divergent parts. 



Fo = 4 log d. 



(4.21) 



The last property of the expression (4.13), which we want to notice here, is its symmetry under inter- 
changing X and y, or equivalcntly. Mi and Jf. 



{ai ...Us] Pi. ..Pi\a,l3,j)g = {-!)'' (f3i . . .^;;ai . ..as\l3,a,j) 

This relation is equivalent to the symmetrization — > — e^ in the formula (3.23) 
3. In the d.s.l. e — s- we may put 

y = -^^ x--^ + V2uo + ..., 



a- 



(4.22) 



(4.23) 



and the equation (4.11) arises. The scaling behaviours of the momenta M^, and d are 

j,__2-(3'=/2+i)e(3fc+i)/2j^ + 4^^ 
Mk -> -2('=-i)/2£-('=-i)/2((2fc - 1)!!4. - 4i), 

d _ 23/2£-3/2 



(4.24) 



Thus, only terms of the highest order in ai that are independent on Jk survive in the d.s.l. when 
the expression ( 4.13| ) converts into the answer for the Kontsevich model (4.£). Then the coefficients 
(ai . ■ .as] {nothing} |a, 0, 7)^ coincide with the Kontsevich intersection indices (tq,j . . . Ta^)g. In | p2[ an 
iterative procedure was proposed for finding coefficients of the expansion ( 4.13| ) ; all these coefficients were 
found in the genus 2 (for g = 0, 1 sec (l8|). It was proved that coefficients of the highest order in 
coincide in a proper normalization with the Kontsevich indices. 



4.2 Relation between momenta and d.m.s. variables. 



We are going to express (4.13) in terms of the quantities standing in the R.H.S. of ( 3.23| ). 

At first, let us expand both momenta Mk, Jk and the restriction equations ( 4.18| , 4.19 ) in terms of 
A-variables, where 77 = ^/a{ e + e 



Then, for the endpoints of the cut, we have: 

x==2Va + ^, y = -2y/a + l3, 



(4.25) 



where ^ and P themselves are some polynomials in the higher momenta Mi and Ji with i,j > 0. Thus, 
after a little algebra we shall obtain, say, for the moment M^: 



Mfc = — tr 
N 



(c^) 



1)2 _ 



fe+1/2 



(e^ + l)2 - ^C^ 



1/2 



(4.26) 
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(For Jfc the expression is just the same with interchanging the powers k + 1/2 and 1/2 for the two terms 
in the denominator.) 

It is convenient now to introduce new momenta: 

Mfe : 



\r. 



N (?7 - x)'=+i/2 ' 
J, = 4tr ^ 



N {r^-y)k+i/2' 
that are related to the initial ones via the following relations: 

Mfc = Mfc_i + 4,2 + d{Mk + Sk,i), 

Mo = Jo = (Mo - Jo)/d. 
Then for these new Alk we have 



1 1 (e^ + l)e^'^ 
■ tr ■ 



1 







1/2 



TV " ,7^'^- (e^ - 



1 - 



\/^(o^-l)2 



fc+1/2 ■ 



The expansion in ( 4.29 ) goes over the terms 



1 (e-^ + l)e'^^ 
Hab — T7 tr ■ 



N (e^-l)2a+i (e^ + l)2''' 

where & > 0, a > fc. 

Let us prove now that Hab can be presented as a linear sum of 

1 



L2a 
R2b 



1 d^" 
tr 



N dX^" e^ - 1' 



1 

tr ■ 



1 



N dX^b c^ + 1' 
i.e., the sum goes only over even powers of derivatives in A: 

a b-l 
Hab^Y.'^abL2^+J2f3i,R2J■ 

This assertion follows directly from the symmetry properties of Hab- 

Hab( — X) ~ —Hab{X), 

= (-l)'+iL,(A)-4o, 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



Thus, Lemma 1.2a is proved. 

Keeping only terms of zero and first orders in traces of A in the expressions for momenta we get: 



Mk 

Jk 

d 



1 1 



, A(fc+1) 



■tr 



■fc+i^V (e^ - l)2fe+i(c^ + 1) 



1 1 e^C^+i 



.fc+iAr*'^(eA_i)(eA + i)2fc+i + ^^i, 
2 



/a {A • — tr 



a N (e-^ - l)(e-^ + 1) 



(4.34) 
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The terms surviving in the d.s.l. are just the ones arising from the term without reductions on the L.H.S. 
of ( 3.23| ). Therefore, we eventually prove (3.5): 



((Tdi • • • TdJ)g = {Td^ ■ . . Tdjg for di + . . . + (i„ = 3.9 - 3 + ; 



(4.35) 



Note that La and Ra are just analogues of the Kontsevich's times r„ = (2n — 1)!!A^"+^. They can be 
transformed into T„ • 2"(n — 1)! in the d.s.l. and in both cases there is no dependence on odd derivatives 
in A. As we have mentioned in Introduction there are two possible ways to do d.s.l. in this model. In the 
first scenario ( 4.23| -4.24) only one set of times {Ln} survives and it turns into the set {T„} in the limit 
£ ^ 0. But if we choose A to be symmetrical, A = diag{Ai, . . . , Xn/2i ^-^i: ■ • • j ^Xn/2}j then another 
limiting procedure is possible where Li = Ri, and each of these sets generates {T„} thus producing a 



square of the integral ( 2.16 ). 

The last note on the reduction procedure concerns the sum over multicomponent reductions on the 
L.H.S. of (3.23). Using the matrix model technique we have an opportunity to distinguish between 
different types of reductions mostly due to the remarkable fact that symmetrization e ^ — ^ — e ^ goes in 
each component separately. Only this property makes A-dependent terms different for, say, {ti{tq)^)o ■ 
(t2To)i and (t2Ti(to)^)o (see Fig. 6) — both these terms appear in the reduction procedure of the genus 
two surface with two punctures. But one of them is due to the two-component reduction and another 
is of one-component type. Evidently, while fixing the number of punctures, n, only terms containing 



products of exactly n traces of A contribute to the L.H.S. of (3.23). 




5 The moduli space A^ij. 



Let us consider an example of modular space Mi^i, i.e. the torus with one puncture. One can immediately 
imagine the copy of the modular figure in Teichmiiller upper half-plane — a strip from —1/2 to 1/2 along 
the imaginary axis bounded from below by a segment of a semicircle of the radius 1 with the origin at 
zero point. In order to get the modular space itself we should identify both sides of the strip as well 
as two halves of this segment being correspondingly on the left and on the right of the imaginary axis 
3?z ~ 0. There are three points where the metric on the moduli space is not conformally flat, namely, 
z ~ i (square point), z = q^'^I'^ (or, the same, e^*'^/'^) (triple point), and z = ioo (infinity point). All 
these points also have a property that each of them is stable under the action of some operator from the 
modular transformation group. For triple point the subgroup of such operators has the third order, for 
square point it is of order 2, and for the infinity point - of an infinite order. 

Thus, the modular space M.\^\ is an orbifold of the (open) upper half-plane. It was namely this way 
how Harer and Zagier j2^ introduced virtual Euler characteristics for such spaces. And it was Penner 
p3| who found a simple one-matrix hermitian model with the potential log(l + X) — X which generated 
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these characteristics. In the Penner approach a factor 1 is assigned to each edge (instead of an 

arbitrary length in the Kontsevich case). Then for an arbitrary Mg,n there is one-to-one correspondence 
between the cells of the simplicial decomposition of the open moduli space Mg^n and the graphs of the 
Penner model. Symmetrical coefRcients for cells and corresponding graphs coincide. Then the virtual 
Euler characteristic be calculated using the formula: 



^ #AutG" 

cells " 
(Graphs) 



(5.36) 



where no is the codimension of the cell in the simplicial complex. 

In the case of Mi^i the triple point graph corresponds to the higher dimensional cell, the square 
point graph - to the cell of codimension 1. In the complex there is also an infinity point of the lowest 
dimension, but there is no graph corresponding to it. Thus, for the virtual Euler characteristic we get 

'^M = \ ■ (-1)° + ^ • (-1)' + - • (-1)' = (5-37) 
3 2 oo 6 

Let us now consider the same case, but already in the Kontsevich-Strebel parameterization. We know 
that there are three types of diagrams depicted in Fig.4a-c. The case of Fig. 4a where all k are greater 
than zero corresponds to the cell of the higher dimension. In A4i^i it is a domain where = -P/^j 

i.e., it is an interior of the equilateral triangle. Note that due to two possible choices of orientation there 
are two such congruent cells. The next case is when one of li is equal zero (Fig. 4b). Taking various U to 
be zero we drive to the boundary of the previous case, i.e. , we get open intervals lying on the boundary 
of the triangles. But it is not the whole boundary as yet — there remains one point at the summit of 
the triangles and it corresponds to the last case. Fig. 4c, where two of li are equal zero, that is a point 
of reduction. The unique reduction of the torus with one puncture is the sphere with three punctures 
whose modular space consists from only one point. 

Now let us draw this simplicial complex 7^ i graphically (Fig. 7). We are to identify the opposite 
edges of the parallelogram thus obtaining the torus. This torus complex consists from two open triangles 
(Fig. 4a), three edges separating these triangles (Fig. 4b) and the unique point of reduction - the vertex 
(Fig. 4c). The centers of the triangles marked by small discs correspond to two copies of the triple point 
and centers of edges - to three copies of the square point (small circles) . We have six copies of the original 
modular figure on this torus, one of them is hatched. 




Fig. 7. A simplicial complex for Mi.i in the Kontsevich picture. 

Thus we result in the conclusion that in the Kontscvich's parameterization the modular space Aii^i 
is the orbifold of a torus with parameters (1, e^^^^) which possesses an internal symmetry group Gi_i 
of the sixth order: 

Afi4=TVGi,i. (5.38) 

This torus is a totally flat compact space. And here is a point which is different from the Penner 
construction of orbifolds of the upper half-plane, because there all infinity points are of infinite order. 
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and here the order of this point is obviously finite! It means that for this case formula (5.37) will change, 
and using ( ^.36 ) we should add 1/6 to (5.37) thus obtaining zero for our new "virtual Euler characteristic" 
in the Kontsevich picture. 

To complete this geometric part, note that we can think about the torus as a fundamental domain 
of the subgroup T2 of the modular group. This domain is depicted on Fig. 8 and it again contains six 
copies of the modular figure. Black discs mark the positions of triple points and small circles - the ones 
of square points. If we identify the left half-line with the right half-circle and vice versa we shall obtain 
the torus (if we do not care about conformal properties of this transformation at the infinity point) . 




Fig. 8. The fundamental domain for subgroup T2 of the modular group. 



Let us turn now to our basic formula ( 3.23| ). First, using diagram technique for the matrix model 



(l.C) it is easy to get the answer (after substitution A = e ). Combining all terms we obtain: 

, 3e2^-l 



6(e2^ - 1)3' 



(5.39) 



and we need to express it in terms of derivatives (4.31). Note that formula (5.39) can be obtained from 



the expansions ( 4.20 ) and ( 4.34 ) substituting a —a/2. After a little algebra we get an answer: 



48a dX^ 



+ 1 



1 

12^ 



(5.40) 



The first term gives us the proper value of {ti)i = 1/24. As for the second term, it originated from the 
"sum over reductions" and the only reduction of the torus is the sphere with three punctures, for which 
(tq)q = 1. We see that the sum over reductions gives an additional fractional factor 1/6, but now we 
know the nature of it. In the simplicial complex (Fig. 7) there are six copies of the modular space Mi^i 
and only one of the infinity point. So we see that we just have "one sixth" of this point contributing to 
the expression (3.23) in this order in g and n. 

Now we are able to understand the structure of d.m.s. for Mi.i- In the Kontsevich parameterization 



we use the form fl (2.6) in order to evaluate the volume of the corresponding modular space. Since the 
intersection indices coincide for both continuum and discrete cases, it does not matter how we calculate 
the total volume of the torus : either by standard continuum integration or by doing a sum over points 
of integer lattice, each taken with unit measure. For the torus with the perimeter equal p (which is always 
even) there are exactly (p/2)2 points from d.m.s. lying in Ti^i = T^. Thus the total volume per one copy 
of the initial moduli space is (p/2)2 divided by the number of copies, i.e., /2A in our case. 



6 Conclusions 

Let us summarize the obtained results and discuss some unresolved problems. 

1. We hope that the established correspondence between the model (l.C) and the discretized moduli 
spaces (Theorem 1.2) may be useful for the understanding of the structure of A^g,™- Here we can select 
the following topics: 
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First, in the Kontsevich-Strebel parameterization it seems that the compactification of the moduh 
space Mg,n is not by Dehgne-Mumford, since (as the example of Mis demonstrates) all points of 
singular curves have symmetry groups of finite orders. In the standard Teichmiiller picture all such 
points have infinite order symmetry groups thus giving zero contribution to the corresponding virtual 
Euler characteristics. It seems also true that in the Kontsevich-Strebel picture each moduli space Mg^n 
possesses a covering manifold, Tg^„, i.e., Mg.n = Tg ^/^g, where Tg is a symmetry group of a finite order. 

Second, due to a possible nonzero curvature of the covering manifold intersection indices ((. . .))g may 
differ from the original (. . .)g for lower orders '^di < d. Therefore, these curvature points (submanifolds) 
could be singular points for the Poisson structures on Mg^n- Here the problem of extracting of symplectic 



leaves appears |26 , |27|. This problem is closely related to the problem of quantization of these structures; 
an attempt in this direction was made by the author , where a quantum group structure was proposed 
in case of There the exceptional representations played an important role. However, all these 

questions are still on the stage of formulation rather than resolving. 

2. A remarkable but a little bit mysterious relation (1.3) establishes a direct bridge from the model 
(l.C) and, therefore, a one-matrix model, to the Kontsevich matrix model. We stress that there is no 
limiting procedure and one can in principle invert the relations (1.3) and (l.A). Moreover, the intertwining 
operator A (1-3) is similar to the free field representation operator in the conformal field theory and both 
one-matrix model and the Kontsevich model are obviously interacting models. At present we are unable 
to give any reasonable interpretation for the very existence of such relation. 
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Appendix A Proof of Theorem 1.3. 

In this Appendix we find and solve constraints on the partition function (l.C) in terms of times corre- 
sponding to d.m.s. 

A.l Algebra of times. 

We express partition function of (l.C) in times 

if = — tr 7 -7 — r—x • (A.l) 

We shall find "fusion rules" for times (A): 

1 r)*^ 1 



fc! 9A/ e^^ ± 1' 



An expansion forniulci for — rrrr 



1 _ 1 _ 1 , B2,n+2 ,2m+l q\ 

"A 2 + ^ (2m + 2)!^ ^^-^^ 



29 



where Bm are Bernoulli numbers. Expansion of -^rrfi is quite the same in the vicinity of the pole A = in, 
but we also need its expansion in A around origin. Using = 77x37 ~ ^2a_x have 

m— 

Taking derivatives for both expressions we get that all odd times are strictly symmetrical under changing 
the sign of A: A ^ —A: 

We find "merging relations" for tJ(A) observing that "negative" times contain only pure poles of 
k + 1-th order in A. From here we can deduce for odd times: 

j.± (y\,± (\\— + i± (\\^\^ B2(n+m-k+l) ^2fc+l(^) 

Hn+A^}H^+iK^) - ± '2(n+™)+3^^^^Z^2(n + m-A; + l) (2m-2fc)!(2n+l)! 



fe=0 



^2(n + m-p+l) (2n-2p)!(2m + l)!' ^ 



and for mixing relation: 



t2n+i{^}Hm+A^) - A.2(n + m-fc + l) (2n-2fc)!(2m+l)! ^z^+i^-^^ 



>p (2^(»+"^-'-+^) - 1) 

^2(n + m-fc + l) (2m-2fc)!(2n+l)! ^z^+i^ ^ ^ ^ 
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A. 2 Constraints on partition function (l.C). 

In this chapter we derive an algebra of constraints imposed on the partition function of the model (l.C) 
in terms of times t^f.. We start with the matrix integral over hermitian N x N matrices X: 



w{e^) = log 



' / DX exp -aN tr {\AXAX + \ log(l - X) + \X) ' 



A=e^. (A.7) 



jAat}. Changing 



/ DX exp-aiVtr {\kXAX~\X'^) 

Without loss of generality we may suppose matrix A to be diagonal, A = diag { Ai 
the variable X ^ -K-^l'^Xk-^l^ - 1 we reduce the upper integral in ( [A.7| ) to a standard integral with 
an external field of GKM type. The integral in the denominator of ( |A.'^ can be easily done, as a result 
we have 



LI (A) = log 



(-l)-^(detA) 



2 e 



■trA^ 



n5.,(A.A,-l)-V2(-)-V2 
exp-aiVtr QaT^ + i log AT + i(A + A"i)X 



(A., 



Doing Itzykson-Zuber integration we get rid of the angular variables, and only eigenvalue integration 
remains (we omit irrelevant numerical factors): 



log 



N 



(A.9) 
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Let rjij = {e^ + e Then we can derive explicitly the Schwinger-Dyson (SD) equations for 

^(e^) = J exp|-^tr[xV2 + logX + X(e^+ e-^)]| (A.IO) 

in terms of 77^. Let (•) mean averaging with the exponential measure taken from (A.IO). We have 
J [I?X]^exp|-^tr [xV2 + logX + X(e^+ e-^)]| 



aN 



{x,j + [X %J +7]ij) = 0. 



(A.ll) 



Taking into account that (xij) = ^-^-^—^i^^) and doing one additional external derivative in we 
obtain SD equation for J^(A) (Q): ' 

1 d d 



(aiV/2)2 diijk 



6^k I 1 
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T{e^) = 0. 



(A.12) 



Using the method described in jl^ we reduce ( A.12| ) to the equation in terms of eigenvalues -qi of -q that 



are equal to e ' 
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or, equivalently, 
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We are interested in the set of equations for e™''^' related to 

N 



Commuting these extra factors with the differentials in ( A.13 ), we eventually get: 
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First, we know that there are no poles in u'(A) of the form x^jX^ — since the original expression 



( [A.7| ) is nonsingular at these points. It means that all such terms should factorize into finite sums of 
times and iJ(A). So far, we deal first with the part of (A.15) which does not contain derivatives in 
A. Tedious algebra demonstrates gentle cancellations of all unwanted terms and gives as a result a very 
simple answer: 
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Next, we already prove from other viewpoint in Chapter 4 (and are able to prove the same from explicit 
analysis of times dependence) that w{X) depends only on even times ij^., k > 0. On an intermediate 
stage we get: 
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(A.17) 



From this expression we already can select a part standing by some chosen i^fe_|_i(Aj). Every such a part 
generates some linearly (but not algebraically) independent constraint on e™'^). The only trouble is with 
the middle term with one derivative originated from the "integral" term in SD equations. We treat it 
now in details, since it is the only one which needs some trick to deal with. 



We start with an identity 
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which is actually due to the symmetry of the expression in x and y both under transformation x —x 



and y — > —y. For the L.H.S. of ( A. 18 ) we have: 
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where d^ „^ € C are some coefficients. We multiply the R.H.S. of ( A.l^ ) and the L.H.S. of ( A.19| ) by 
e^ — e^ and e^ — e~^, respectively, in order to eliminate the prcfactors in front of the derivative terms, 
and then subtract one expression from the other. It gives 
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From this relation it is already easy to find that 
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Combining all terms from ( A.17 ), we shall obtain a set of conditions on e^'^'^^ of the form: 

oo cso 



fe=0 



fe=0 



which are valid for any j. We assume that all traces of matrix A are independent. Note, however, that 
from this it does not follow that all times and t^j,^]^ are linearly independent. There exists a formula 
that connects these two sets of times by re-expansion of positive times via negative ones using the Taylor 
expansion and the shift relation tf{X) = —tj{X + in). But as far as we are looking for an expansion of 
w{X) over additional parameters N and a, and in each fixed order in TV and a the time dependence is 
polynomial, then it follows that this expansion is unique (the connection formulas between positive and 
negative times are obviously non-polynomial). So, we treat all operators as independent generators 

of the constraint algebra for e As usual, the first generator, Lf, is somehow selected from the whole 
set and we give it separately: 
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and for s > 
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Analogous formulas for -£/2s+i can be obtained changing all f Jj, — and 

A. 3 Algebra of constraints Lf^^^. 
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We derive now commutation relations for generators ^25+1- Tedious but again direct calculations show 



(A.24) 



i.e. in spite of the fact that generators and L. contain derivatives in both positive and negative times, 
two halves of the algebra factorize. 

Let us consider the algebra of i^. We have 
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After an upper triangular transformation of generators: 
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we arrive to the standard Virasoro algebra : 
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We are interested in the time transformation corresponding to the change of generators (A.26). In 
fact, by analysis in orders of a and N, we know that in order by order calculations the lowest term is 
the one from ( |A.23| ) that is equal to So we look for such deformed times t, in which this term 

preserves its form in Lg. It means that 
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A solution to this set of equations is provided by lower triangular transformed times: 
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(Note that most of these transformations are based on an identity: 
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which also generates famous relations for Bernoulli numbers.) 

In fact, there is no great simplification of the formulas for Lg in terms of these new times. We present 
here only the expression for the part of L^, s > 0, which is linear in derivatives: 
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Eventually, ij^^j^ were combined with times t^s_(_i(A). in (A) standing with Ls are 
tn{\) = <^„+i(A) - E 1^ 7^4„-2fe-i(A), n > 0, 
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2fc + 2 (2fc)! 



or, in terms of the Bernoulli polynomials: 
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where B„(a;) = ^^^^ (^) Bo = I, Bi = -1/2, B2 = 1/6, B4 = -1/30, Be = 1/42, = -1/30 

Bio = 5/66, etc, S2n+i = for n > 0. 

We also present two formulas showing how "non-tilde" times are expressed via i„: 
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Let us now turn again to the "old" (non-tilde) times. We use expressions ( |A.5| ) and ( [A. 61 ) in order 
to simplify expression (A. 22). For example, we consider the term with two "minus" derivatives which 
occupies third to fifth lines in (A. 22). It originates ft'om the term proportional to t(j^(A) in the expansion of 
^2n+iWt2m+iWto W ■ ^^^t expand t2„i+i{^)to W and after that merge it with ^2ti+i('^)i then it 

appears that it simplifies drastically and converts into 
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It is also worth to note that after shift (A. 26) most of "tails" in formulas (A. 22) and (A. 23) disappeared 
and as a result we have simplified expressions for 
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and for L,, s > 0: 
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2(ji+s) 

We checked explicitly that these generators do satisfy two halves of Virasoro algebras Vir+ : 
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where Ls = and are obtained from Lg by the interchange tj, ^ ~^2s ^^^^ ^ 
It is worth noting that expression ( A.33D can be rewritten in the form 
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where ah derivatives are assumed to act on the right (that ehminates the constant term). This formula 
gives us a hint how to simphfy further the expressions ( |A.33 ) and (A.34). 

Let us look for such a canonical transformation of times and their derivatives that shifts times by 
some linear combination of derivatives plus constant terms and leaves time derivatives unchanged: 
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This immediately gives us 
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Making now the last substitution: 
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we eventually obtain a simple expression for the Virasoro generators (A. 33) and (A.34) in terms of T*: 
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But these Virasoro conditions are nothing but the Virasoro constraints of the Kontsevich matrix model! 
(see Q). Their solution is well-known. It satisfies the KdV equations and was elaborated to the third 
and the forth orders in Therefore, there exists a canonical transformation of variables of the initial 
matrix model (A_/7) that reduces it to two copies of the Kontsevich model. Taking into account that the 
vacuum is invariant under these transformations we obtain for (A. 7): 
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where A is found in ( A.39| ), and C(aN) is a function independent of times, such that both sides of ( |A.43 ) 
are units when = 0, 

nttn) ^ mn+l) ^ ^ , (A.44) 

where ^2n+i are odd times of the Kontsevich model and F(^2n+i) is just the partition sum of the Kont- 
sevich model. Therefore we have proved the main assertion of the Theorem 1.3. 

A. 4 Perturbative solution for {{rd^ . . -Td,^))- 
Let us consider an expansion of w{X) in a and N: 
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In terms of times the expansion coefficients are (aN)'^ " and we have an asymptotic expansion 
of the form: 
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Here, taking the times ^2^+1 = (^^ + ^)-T^ in order to compare with the answer for the Kontsevich 
model (0), we have 
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Appendix B The explicit solution to A4 



2,1 



Here we shall present the form of formula ( |3.23 ) for the case of genus two moduli space with one puncture. 



In paper |12| the explicit form of genus two partition function in terms of momenta was found: 

181 181 5 I8IJ2 I8IM2 



F2 = 



480Jld^ 480Af2d4 iQj^Mid'^ A&OJfd^ 480Affd3 
3J2 3M2 11J| 11M| 



+ 



J2M2 43Af3 43J3 21J| 21M| 



GAJfM^d^ 192Mfd^ 192j3d2 160Jfd 160Mfd 
29J2J3 29M2M3 35J4 35M4 

128 Jfd ^ 128Mfd ^ 38Ufd ~ 384Mfd' ^ ' -* 
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In order to investigate the modular space M2,i it is enough to use the expansions (4.34), because we 
must keep only terms of the first order in traces. Then the only thing we need more is to express the 
quantities 



,A(fe+l) 



Pk 

qk 



l)2'=+i(e^ 

gA(fc+l) 



(e^ - l)(e^ + 1) 



(B.2) 



via the derivatives La and Ra (4.31). We omit all lengthy calculations and present here only the final 
answer. After replacing a —a/2 we remain with 
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